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Lattice gluodynamics at negative g 2 



O 

o 

(N 

> 
O 

> 
On 
(N 
O 
O 

o , 

> . 

I 

^ ■ 



X 



L. Li and Y. MeuriceEJ 

Department of Physics and Astronomy 
The University of Iowa 
Iowa City, Iowa 52242 
USA 

(Dated: February 1, 2008) 

We consider Wilson's SU(N) lattice gauge theory (without fermions) at negative values of j3 = 
2N/g 2 and for N—2 or 3. We show that in the limit j3 — > — oo, the path integral is dominated 
by configurations where links variables are set to a nontrivial element of the center on selected 
non intersecting lines. For N = 2, these configurations can be characterized by a unique gauge 
invariant set of variables, while for N = 3 a multiplicity growing with the volume as the number 
of configurations of an Ising model is observed. In general, there is a discontinuity in the average 
plaquette when g 2 changes its sign which prevents us from having a convergent series in g 2 for this 
quantity. For N = 2, a change of variables relates the gauge invariant observables at positive and 
negative values of j3. For N = 3, we derive an identity relating the observables at f3 with those at f3 
rotated by ±27r/3 in the complex plane and show numerical evidence for a Ising like first order phase 
transition near /3 = —22. We discuss the possibility of having lines of first order phase transitions 
ending at a second order phase transition in an extended bare parameter space. 

PACS numbers: 11.15.-q, 11. 15. Ha, 11.15.Me, 12.38.Cy 



I. INTRODUCTION 

It has been known from the early days of QED that 
perturbative series have a zero radius of convergence 0. 
This has not prevented Feynman diagrams to become an 
essential tool in particle physics. However, perturbative 
series need to be used with caution. The divergent nature 
of QED series was foreseen by Dyson as a consequence 
of the apparently pathological nature of the ground 
state in a fictitious world with negative e 2 . Like charges 
then attract and pair creation can be invoked to produce 
states where electrons are brought together in a given 
region and positrons in another. Dyson concludes that 
as this process sees no end, no stable vacuum can exist. 

For Euclidean lattice models, related situations are en- 
countered. For scalar field theory with \<f> interactions, 
configurations with large field values make the path inte- 
gral ill-defined when A < (provided that no higher even 
powers of (f> appear in the action with a positive sign and 
that the path of integration is not modified). Modified 
series with a finite radius of convergence can be obtained 
by introducing a large field cutoff 00- We are then con- 
sidering a slightly different problem. In simple situations 
, it is possible to determine an optimal value of the field 
cutoff that, at a given order in perturbation, minimizes or 
eliminates the discrepancy. For non-abelian gauge theo- 
ries in the continuum Hamiltonian formulation, the sub- 
stitution g — > ig makes the quartic part unbounded from 
below and the cubic part non-hermitian. 

It should be noted that in quantum mechanics 0, it 
is possible to change the boundary conditions of the 
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Schrodinger equation in such a way that the spectrum 
of an harmonic oscillator with a perturbation of the form 
ix 3, or —x 4 stays real and positive. The procedure can 
be extended to scalar field theory in order to define a 
sensible icf) 3 theory Q . Even though conventional Monte 
Carlo calculations would fail for these models, complex 
Langevin methods can be used to calculate Green's func- 
tions 0- 

In the case of lattice gauge theory with compact gauge 
groups, the action per unit of volume is bounded from 
below and there is no large field problem. Consequently, 
these models have well-defined expectation values when 
g 2 < and we can consider the limit g 2 — > CP. In this ar- 
ticle, we discuss the behavior of Wilson loops for SU(N) 
lattice gauge theory with g 2 < 0. This work is motivated 
by the need to understand the unexpected behavior of the 
lattice pertubative series for the I x 1 plaquette calculated 
up to order 10 [3, la DU • An analysis of the successive 
ratios 0, may suggest that the series has a finite 
radius of convergence and a non-analytic behavior near 
P ~ 5.7 in contradiction with the general expectations 
that the series should be asymptotic and the transition 
from weak to strong coupling smooth. 

We consider here pure (no fermions) gauge models with 
a minimal lattice action [13| . For definiteness this model 
and our notations are defined in section [H] The extrema 
of the action are discussed in II I II and enumerated for 
5(7(2) and SU(3). We then discuss (section ITvT) the case 
of SU(2) and show that planar Wilson loops with an area 
A (in plaquette units) pick up a factor (— 1) A when g 2 
becomes negative and the behavior for g < is com- 
pletely determined by the behavior with g 2 > 0. As the 
Wilson loops are nonzero when g — > + , the ones with 
an odd area have a discontinuity which invalidates the 
idea of a regular perturbative series. 

The case of SU(3) is discussed in section Ivl where we 
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derive identities involving Wilson loops calculated with 
a coupling rotated by ±27r/3 in the complex plane. Con- 
sequently, for TV = 3, we have no a-priori knowledge re- 
garding the behavior of Wilson loops when g 2 < 0. We 
report numerical evidence for a first order phase transi- 
tion near f3 = 6/g 2 — — 22 using methods similar to Ref. 
|l4j . The implications of our findings are summarized in 
the conclusions. 



II. THE MODEL, NOTATIONS 

In the following, we consider the minimal, unimproved, 
lattice gauge model originally proposed by K. Wilson 
|13| . Our conventions and notations are introduced in 
this section for definiteness. We consider a cubic lattice 
in D dimensions. A SU(N) group element is attached 
to each link I and U denotes its fundamental represen- 
tation. Up denotes the conventional product of Ui (or 
hcrmitian conjugate) along the sides of a 1 x 1 plaquette 
p. The minimal action reads 



S = f3j2(l-(l/N)ReTr(U p j) 



(1) 



with f3 — 2N/g 2 . The lattice functional integral or par- 
tition function is 



n / du i e ~ 



(2) 



with dUi the SU(N) invariant Haar measure for the 
group element associated with the link I. The average 
value of any quantity £3 is defined as usual by inserting 
13 in the integral and dividing by Z. 

In the following, we consider symmetric finite lattice 
with L D sites and periodic boundary conditions. For 
reasons that will become clear in the next sections, L will 
always be even. The total number of 1 x 1 plaquettes is 
denoted 



= L D D(D - l)/2 



Using 



f = -(l/X p )lnZ 
we define the average density 



(3) 



(4) 



P(/3) = df/d/3 

= (l/tf p ) (^(1 - (l/N)ReTr(U p ))^ . (5) 

In statistical mechanics, / would be the free energy den- 
sity multiplied by (3 and P the energy density. In analogy 
we can also define the constant volume specific heat per 
plaquette 



III. THE LIMIT f3 ->• -oo 

In the limit — » — oo, we expect the functional in- 
tegral to be dominated by configurations which maxi- 
mize YlpO- ~ (l/N)ReTr(Up)). In the opposite limit 
(/3 — > +oo), the same quantity needs to be minimized 
which can be accomplished by taking Ui as the identity 
everywhere. 

We first consider the question of finding the extrema of 
TrU . For our study of the behavior when /3 — > — oo, we 
are particularly interested in finding absolute minima of 
TrU. Using TrU = Tr(VUV^) for V unitary, U = e lH 
with H traceless and hermitian, and Ve lH V^~c lVHvt , 
we can diagonalizc H and write 



JV-l 



N-l 



ReTrU = cos(^) + cos(^ fc) 



The extremum condition then reads 



JV-l 



sin(^i) + sin( ^ fc) = 



(7) 



(8) 



CV = -f3 2 dP/dp 



(6) 



for i = 1, . . . iV — 1. The trivial solution is all 4>i = 0. We 
then have ReTrU = N which is an absolute maximum. 

For N = 2, we have only one nontrivial solution <f>\ — 
7r, which corresponds to the nontrivial element of the 
center U = — 1. We then have ReTrU = —2 which is an 
absolute minimum. 

For N — 3, we have 5 nontrivial solutions. Two cor- 
respond to the nontrivial elements of the center (0i = 
<f>2 = ±27r/3). The matrix of second derivatives has two 
positive eigenvalues and these two solutions correspond 
to a minimum. We use the notation f2 = e l27T ^ 3 l on the 
diagonal. We have ReTril — —3/2, which we will see 
is an absolute minimum. The other three solutions are 
(4>i = 7r, 02 = 0), (0i = 0, cfe = 7r) and (0i = n, (j> 2 = n). 
They correspond to elements conjugated to diagonal ma- 
trices belonging to the three canonical SU(2) subgroups 
with the SU(2) element being the non trivial center el- 
ement. These three solutions have matrices of second 
derivatives with eigenvalues of opposite signs and corre- 
spond to saddle points rather than minimum or maxi- 
mum. In the three cases ReTrU = — 1. 

For general N, it is clear that we can always find at 
least one group element U such that ReTrU is an ab- 
solute minimum. In particular, for N even, U = — 1 
is such a group element, with ReTrU = —N (the indi- 
vidual matrix elements must have a complex norm less 
then one). For N > 3 and odd, it is easy to check that all 
<f>i — [N — l)n/N is a solution of the extremum condition 
Eq. p. For this choice, ReTrU = -N\ cos{{N-1)tt/N)\ 
which is clearly negative and tends to —N as N becomes 
large. This solution (the element of the center the closest 
to —1) gives an absolute minimum of TrU for N = 3 and 
we conjecture that it is also the case for larger N. 

We can now obtain an absolute minimum of the action 
if we can build a configuration such that ReTrU p takes 



3 



its absolute minimum value for every plaquette. This 
can be accomplished by the following construction. In 
the Appendix, we show that (at least for for D < 4) and 
for L even, it is possible to construct a set of lines on 
the lattice such that every plaquette shares one and only 
one link with this set of lines. We call such a set of links 
£. One can then put an element which gives an absolute 
minimum of ReTrU on the links of £ and the identity on 
all the other links. For SU(2), there is only one possible 
choice that minimizes ReTrU, namely U = — 1. For 
SU(3), there are two possible choices U = £1 or U = W. 
We emphasize that the construction only works for L 
even. If L is odd, there will be lines of frustration in 
every plane. 

The set of links £ is not unique. Starting with a given 
set, we can generate another one by translating the lines 
by one lattice spacing or rotating them by ir/2 about the 
lattice axes. By direct inspection, it is easy to show that 
for D = 2 there are 4 such a sets of lines while for D = 3 
there are 8 of them. 

Enumerating all the gauge inequivalent minima of the 
action at negative (3 for arbitrary D and N appears as 
a nontrivial problem. In the rest of this section, we spe- 
cialize the discussion to N = 2 or 3. In order to dis- 
criminate among gauge inequivalent configurations, it is 
useful to make the following (gauge invariant) argument: 
in order to have an absolute minimum of the action, for 
every lxl plaquette p, the product U p (n) of the Ui 
along p starting at any site n of p, is a nontrivial ele- 
ment of the center. Under a local gauge transformation, 
U p {n) -> V(n)l7p(n)V(n)t = U p (n) since U p (n) com- 
mutes with any SU (N) matrix. For the same reason, 
changing n along the plaquette amounts to a VUV^ con- 
jugation and has no effect on the center. Consequently, 
configurations with a different set of it = {U p } are not 
gauge equivalent. One can think of it as a set of electric 
and magnetic field configurations. 

For 5/7(2), there is only one, uniform, set it where all 
the elements U p = — 1. For D = 2, this can be realized 
in 4 different ways by putting — 1 on the 4 distinct sets 
£. These 4 configurations are all gauge equivalent. The 
gauge transformations that map these 4 configurations 
into each others can be obtained by taking V = — 1 on 
every other sites of the lines of £. For D = 3, it is also 
possible to show that the 8 configurations that can be 
constructed with a similar procedure can also be shown 
to be gauge equivalent. The gauge transformation can 
be obtained by taking V = — 1 on every other sites of 
the lines of £ pointing in two particular directions and 
in such way that one half of the lines created by the 
gauge transformation associated with one direction "an- 
nihilate" with one half of the lines created by the gauge 
transformation associated with the other direction. Wc 
conjecture that in higher dimensions, the configurations 
that minimize the action for SU(2) are also related by 
gauge transformations. 

For ST/ (3) the situation is quite different because for 
every link of a particular £, we have two possible non- 
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FIG. 1: The average action density P{0) for £77(2). 

trivial element of the center. Since there are "N p = 
L D D(D - 1)/2 plaquettes on a L lattice and one link of 
£ per plaquette, shared by 2(D — 1) plaquette, we have 
DL D /4 links in any £. Picking a particular £, it possible 

to construct 2 DL / 4 distinct it. Consequently, there are 
at least 2 DL / 4 gauge inequivalent minima of the action 
for SU(3) . Note that 2 DL / 4 always is an integer for L 
even, which has been assumed. In the case D = 4, the 
degeneracy is simply 2 which is the same as the number 
of configurations of an Ising model on a L 4 lattice. 

In summary, we predict a discontinuity in P as g 2 
changes sign. In the limit f3 — > +00, we have P — > 0, 
while in the limit /3 — > -co, we expect P — > 2 for N 
even, and 1 + | cos((7V - 1)tt/N)\ for N odd. 

IV. N=2 

In this section we discuss SU (2) gauge theories at neg- 
ative (3. The basic idea is that it is possible to change 
fiReTrUp into —(3ReTrU p by making the change of vari- 
ables Ui — > — Ui for every link / of a particular £. Since 
— 1 is an element of SU(2) and since the Haar measure 
is invariant under left or right multiplication by a group 
element, this does not affect the measure of integration. 
Consequently, we have 

Z(-0) = e 2pys 'Z{P) (9) 

Taking the logarithmic derivative as in Eq. JSJ), we obtain 

P(/3)+P(-/3) = 2. (10) 

This identity can be seen in the symmetry of the curve 
P(J3) shown in Fig. [TJ The validity of Eq. JTDJ can be 
further checked by calculating the difference 

A(/?) = |P(/3)+P(-/?)-2| (11) 

which should be zero except for statistical fluctuations. 
Fig. |2 illustrates this statement and shows that the sta- 
tistical errors of our calculations are of order 10~ 4 or less. 
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FIG. 2: A(/3) defined in Eq. JTTJ versus /3. 

The relation between P(±/3) of Eq. QlOfl together with 
the assumption that P(+oo) = 0, is in agreement with 
the statement made in section ITTTl that P seen as a func- 
tion of g 2 — 2JV//3, jumps discontinuously by 2 as g 2 
becomes negative. This invalidates the idea that P could 
have a regular expansion about g 2 = with a non-zero 
radius of convergence. 

This relation can also be used in the opposite limit 
and expanded about f3 = 0. The odd terms cancel au- 
tomatically. The even terms of order 2 and higher add 
and cannot cancel. Consequently, the even coefficients of 
the strong coupling expansion of P((3) and the odd co- 
efficients of the free ener gy s hould vanish, in agreement 
with explicit calculations [lfij . 

The discontinuity at g 2 = can be extended to Wilson 
loops of odd area (in plaquette units). To see this, let us 
consider a Wilson loop W(C) with C a contour that is 
the boundary of an area made out of A plaquettes. For 
simplicity, let also assume that this area is connected and 
has no self-intersections. Under the change of variables 
Ui — > — Ui for every link I of an arbitrary set £, we have 
W{C) -> (-1) A W{C). This follows from the fact that 
for any line, the parity of the number of links of C shared 
with this line, is the same as the number of plaquettes of 
the area in contact with this lines. Since £ shares a link 
with every plaquette, we obtain the desired result. This 
can be summarized as 

{W{C))_ p = {-l) A {W{C)) p (12) 

We can now try to interpret the change of the Wilson 
loop with the area in a term of a potential. We consider 
a rectangular RxT contour C and write 

W(R, T, 0) ee (W(C))p oc e ~ E{R ^ T . (13) 

From Eq. (|12|) this implies 

E(R,-\/3\) = E(R,\P\) +iwR . (14) 

This property can be related to the fact that the config- 
urations of minimum action are invariant under transla- 
tions by two lattice spacings but not under translations 



by one lattice spacing. This also confirms our expecta- 
tion that the hamiltonian develops a nonhermitian part. 

V. N=3 

For N = 3, — 1 is not a group element and the closest 
thing to the change of variables used for N = 2 that we 
can invent is a multiplication by a nontrivial element of 
the center £1 for the links of a particular set £. We then 
obtain 

Z((/3) = eV-QPXpztf) 



e (/3/3) E p (CR<*TrU p -ReTrU p ) 




In the case N = 2, £ is replaced by -1, (•••)« becomes 
1 and we recover Eq. (JSJ. In the case of SU(3), the 
factor (...)^ prevents us from deriving an exact identity 
analog to Eq. i(TT)|) for SU(2). It is however possible 
to obtain an approximate generalization which is a good 
approximation for small (3. Setting (3 = £*cc, taking the 
logarithmic derivative with respect to x and setting x = 
we obtain 

p(0) = i-e + cp(p) + ow . (i6) 

Taking the real part and using 1 + £ + £ 2 = 0, we obtain 

P(/3)+P(C/3)+P(C 2 /?)=3 + 0(/3 6 ) . (17) 

which can be seen as an approximate SU (3) version of 
Eq. 110|l . The cancellation of the terms of order 1, 2, 
4 and 5 occurs independently of the values of the coef- 
ficients at these orders. The absence of contribution of 
order 3 and the presence of a nonzero contribution at or- 
der 6 comes from the fact ^(J that ln(Z)/ s Jl p has a zero 
(nonzero) contribution at order 4 (7). 

As it does not seem possible to obtain P{(3) for (3 real 
and negative from our knowledge at (3 real and positive, 
we have to resort to a direct numerical approach. The 
results are shown in Fig. [3] A discontinuity near (3 = 
—22 is clearly visible. This indicates a first order phase 
transition. 

The metastable branches have been studied following 
the approach of M. Creutz 0] used to study of a fist or- 
der transition for SU(5). As (3 becomes more and more 
negative, the system becomes more ordered but has also 
higher average energy P, the supercooled/heated termi- 
nology may be confusing and will be avoided. 

We have run Monte Carlo simulations on a 8 4 lattice 
at (3 = —22 with 4 different initial configurations. Our 
best estimate of the critical (3 for this volume is -22.09. 
The first initial configuration was completely ordered (in 
the (3 — > — oo sense, with P=1.5) by putting a nontriv- 
ial element of the center on a set of lines £. As we set 
[3 = — 22, we expect to stay on the upper branch and end 
up with P ~ 1.39 (black dots in Fig. for many itera- 
tions. The second configuration was completely random 
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FIG. 3: The average action density P(/3) for SU(3) 



(empty circles) and stayed on the lower branch when (3 
was set to -22, to end up at P ~ 1.34. The third con- 
figuration (empty squares) was initially random, we then 
temporarily set {3 = —27 letting P go up to 1.38, ex- 
pecting to reach lower metastable branch. When (3 is set 
to -22, P stabilized to the lower value 1.34. Finally, we 
prepared a fourth initial configuration (empty triangles) 
by first setting a nontrivial element of the center on a 
given £ and then temporarily setting (3 = — 17 until P is 
near 1.35 expecting to reach the upper metastable part. 
When (3 is finally set to -22, we reach the upper branch 
value P = 1.39. Fig. 0]is quite similar to Fig. 1 of Rcf. 
[l4| and has the same type of crossings. 
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FIG. 5: Distribution of ImTrU for three values of /3. 
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FIG. 4: P as a function of iterations for the four initial con- 
figurations described in the text. 

We believe that the first order transition observed 
above is similar to the one observed^]} for the D = 4 Z% 
gauge theory. This model is dual to a nearest neighbor 
Ising model. In Fig. we show histograms of the dis- 
tribution of ImU below, near and above the transition. 
ImU allows to separate the two nontrivial elements of the 
center Q and f^. As (3 becomes more negative and goes 
through the transition, a broad distribution around de- 
velops two bumps which keep separating and sharpening 



as one would observe in an Ising model. 

The transition can also be seen as a singularity in the 
specific heat defined in Eq. 10 as shown in Fig. [S] As 
expected the height of the peak increases with the vol- 
ume. The location of the transition sightly moves left as 
the volume increases. 
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FIG. 6: The specific heat CV versus j3 near the first order 
phase transition. 

Finally, we would like to compare the decay of the Wil- 
son loop at negative (3 for SU(2) and SU(3). In Fig. [7| 
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we have plotted the Wilson loop (W(l, R)) for these two 
groups. For SU (2) we observe the same decay as at pos- 
itive /3 but with alternated signs as predicted in section 
ITVI For SU(3), the decay is much faster than at the pos- 
itive value of j3 and show that the sign alternates as long 
as the signal is larger than the statistical fluctuations 
(namely for R < 6). 

The sign alternates at relatively low values of \(3\. This 
agrees with the strong coupling expansion which predicts 
a (— |/9|/18)^ x1 ' behavior for (3 negative and small in 
absolute value. 
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FIG. 7: The Wilson loop (W(1,R)) at (3 
and = ±6 for SU(3). 
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VI. CONCLUSIONS 

We have studied lattice gauge theories at negative l3. 
Wilson loops are well-defined and calculable with the 
Monte Carlo method. However, the limits (3 — ► ±oo of P 
differ and an expansion in g 2 cx 1//3 cannot have a finite 
radius of convergence. This statement has been substan- 
tiated for N = 2 and 3, but from the discussion of section 
IIIII it seems clear that it should extend to general N . 

We found a first order phase transition near (3 = —22. 
At this point, it seems unrelated to the known transition 
near (3 ~ +6 and branch cuts in the complex plane dis- 
cussed by J. Kogut [18| . However, a more complete pic- 
ture may appear if we study P for a larger class of action. 
It is conceivable that by introducing a linear combination 
of terms involving larger contours than the lxl plaque- 
tte, multiplied by one free parameter, we could create 



a line of first order phase transition ending at a second 
order phase transition. Another possibility would be to 
add an adjoint term as in Ref. |19|. Finding new second 
order phase transitions would allow us to define a non- 
trivial continuum limit. This could be of interest in the 
context of cosmology and stellar evolution. 
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APPENDIX: MAXIMAL SETS OF 
NON-INTERSECTING OF LINES ON A CUBIC 
LATTICE 

In this Appendix, we consider a D dimensional cubic 
lattice. We restrict the use of "line" to collections of 
links along the D principal directions of the lattice and 
the use of "plane" to collections of plaquettes along the 
D(D — l)/2 principal orientations. In other words, these 
objects are lines and planes in the usual sense, but we 
exclude some "oblique" sets that can be constructed out 
of the sites. 

We now try to construct a set of lines such that every 
plaquette shares one and only one link with this set. It 
is obvious that these lines cannot intersect, otherwise, 
at the point of intersection and in the plane defined by 
the two lines, we could fit 4 plaquettes, each sharing two 
links with the lines. These lines cannot be obtained from 
each other by a translation of one lattice spacing in one 
single direction, otherwise the set of lines would share 
two opposite links on the plaquettes in between the two 
lines. 

For D = 2, the problem has obvious solutions, we can 
pick for instance a set of vertical lines separated by two 
lattice spacings. Using translation by one lattice spacing 
and rotation by n/2, it is possible to obtain three other 
solutions. For D > 2, it is sufficient to show that for 
every plane (in the restricted sense defined above), we 
have a D = 2 solution. As this restricted set of planes 
contains all the plaquette ounce, we would have then suc- 
ceeded in proving the assertion. If such a solution exists, 
it is invariant by a translation by 2 lattice spacing in any 
direction. Consequently, we only need to prove the ex- 
istence of the set of lines on a 2 D lattice with periodic 
boundary conditions. The full solution is then obtained 
by translation of the 2 D solution. If the lattice is finite, 
this only works if L is even, an assumption we have main- 
tained in this article. 

On a 2 D lattice, the lines (as defined above) are con- 
structed by fixing D — 1 coordinates values to be or 
1 and leaving the remaining coordinate arbitrary. For 
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instance, for D = 3, a line in the 3rd direction com- 
ing out of the origin will be denoted (0, 0, A) where A 
stands for arbitrary and means or 1. In general D, 
there are D2 D ~ 1 such lines. Consequently, there are 
D2 D links, each shared by 2(D — 1) plaquettes. There 
arc thus D2 D 2(D - l)/4 = D(D - l)2 D - 1 plaquettes. A 
set of lines which has exactly one link in common with 
every plaquette, has D(D - 1)2 L) - 1 /(2(D - 1)) = D2 D ~ 2 
links in other words it must contain D2 D ~ 3 lines. For 
D = 2, such a set has only one line and there are four 



possible choices. For D = 3, an example of solution is 
{(A, 0, 0), (0, A, 1), (1, 1, A)}. It is not difficult to show 
that there are 8 distinct solutions of this type. For D = 4, 
a solution consists in 8 lines. An example of solution is 

{ (A, 0, 0, 0), (0, A, 0, 1), (0, 1, A, 0), (0, 0, 1, A), 
(1, 1, 0, A), (1, 0, A, 1), (1, A, 1, 0), (A, 1, 1, 1) } . 
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